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Abstract 

(N 

^-H Particle IDentification (FID) is fundamental to particle physics experi- 

ments. This paper reviews PID strategies and methods used by the large LHC 
' ^ ' experiments, which provide outstanding examples of the state-of-the-art. The 

first part focuses on the general design of these experiments with respect to 
-i. PID and the technologies used. Three PID techniques are discussed in more 

IjJ detail: ionization measurements, time-of-flight measurements and Cheren- 

kov imaging. Four examples of the implementation of these techniques at 
the LHC are given, together with selections of relevant examples from other 
experiments and short overviews on new developments. Finally, the Alpha 



rn 

^ Contents 
O 

^ 1 Introduction 0] 

1.1 PID by difference in interaction |3] 

• '-j 1.1.1 Tracking system ID 

rS 1.1.2 Calorimeters g] 

c3 1.1.3 Muon system |6] 

1.1.4 Other particles |6] 

1.2 PID by mass determination |7] 

2 Overview of the large LHC experiments |9] 

2.1 ATLAS and CMS M 

2.1.1 Requirements [10] 

2.1.2 Setup \n\ 



* Christian.Lippmann @ cern.ch 



1 



2.1.3 Tracking and muon systems [TT] 

2.1.4 Particle identification [13] 

2.2 ALICE m 

2.2.1 Requirements [15] 

2.2.2 Setup [16] 

2.2.3 Particle identification [TT] 

2.3 LHCb \M 

2.3.1 Requirements [18] 

2.3.2 Setup M 

2.3.3 Particle identification [19] 

3 Ionization measurements HO] 

3. 1 Energy loss and ionization [20] 

3.2 Velocity dependence [21] 

3.3 Straggling functions [22] 

3.4 PID using ionization measurements [22] 

3.5 Energy resolution and separation power [24] 

3.6 Errors affecting the resolution [25] 

3.7 ALICE Time Projection Chamber [26] 

3.7.1 Field cage and readout chambers [26] 

3.7.2 Front-end electronics [27] 

3.7.3 PID performance and outlook [28] 

3.8 Silicon detectors at LHC [29] 

3.9 Developments for future TPCs [30] 

4 Time-of-flight |3i 

4. 1 Time resolution and separation power [30] 

4.2 Errors affecting the resolution |3T] 

4.3 Resistive Plate Chambers for time-of-Hight measurements .... [32] 

4.4 Detector physics of Resisitive Plate Chambers [32] 

4.5 ALICE time-of-flight detector [33] 

4.5.1 Time resolution [34] 

4.5.2 PID performance and outlook [36] 

4.6 Other TOF detectors [36] 

4.7 Developments for future TOF systems [38] 

5 Clierenkov imaging |39] 

5. 1 Cherenkov ring imaging |40] 

5.1.1 Radiator E] 

5.1.2 Optics gl] 

5.1.3 Photon detection |4T] 

5.1.4 Pattern recognition |42] 

5.2 Angular resolution and separation power [42] 

5.3 ALICE HMPID detector |43] 



2 



5.5 
5.6 



5.4 



5.3.1 Radiator 

5.3.2 Photon detector 

5.3.3 Performance and outlook 

LHCb RICH 

5.4.1 Radiators 

5.4.2 Photon detectors 

5.4.3 Performance 

Other Cherenkov detectors 

Developments for future RICH detectors . 

5.6.1 Proximity-focusing aerogel RICH 

5.6.2 Time-of-propagation counters . . 



mi 
mi 

M 

m 
m 
m 
m 
m 
m 
m 
m 



6 PID in space 



7 Summary and outlook 



m 



1 Introduction 

In addition to tracking and calorimetry, Particle IDentification (PID) is a crucial 
aspect of most particle physics experiments. The identification of stable particles 
is achieved either by analysing the way they interact, or by determining their mass. 
The difference in interaction is primarily used for lepton and photon identification. 
A "traditional" particle physics experiment already incorporates this method in 
its conceptual design, as it is discussed in Sect. |1.1| In order to unambiguously 
identify hadrons, their charge and mass has to be determinded, which is achieved by 
simultaneous measurements of momentum and velocity. This method is discussed 
in Sect.fOl 

The four large experiments at the CERN Large Hadron Collider (LHC) HI and 
their PID strategies are introduced briefly in Sect. [2] Three methods to determine 
the mass of a charged particle and examples of their implementation at the LHC 
experiments are discussed in Sects. [3] (ionization), [4] (time-of-flight) and [5] (Che- 
renkov radiation imaging). A fourth important method, the detection of transition 
radiation, is not discussed, since it is the topic of a separate paper in this issue ||2l. 
In Sect. [6] finally the AMS 02 experiment is briefly described as an example for 
the usage of PID techniques from accelerator-based particle physics experiments 
in space. 

1.1 PID by difference in interaction 

In a "traditional" particle physics experiment particles are identified (electrons and 
muons, their antiparticles and photons), or at least assigned to families (charged 
or neutral hadrons), by the characteristic signatures they leave in the detector. The 
experiment is divided into a few main components, as shown in Fig. [T] where 
each component tests for a specific set of particle properties. These components 
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Figure 1: Components of a "traditional" particle physics experiment. Each particle 
type has its own signature in the detector. For example, if a particle is detected 
only in the electromagnetic calorimeter, it is fairly certain that it is a photon. 



are stacked in layers and the particles go through the layers sequentially from the 
collision point outwards: first a tracking system, then an electromagnetic (EM) and 
a hadronic calorimeter and a muon system. All layers are embedded in a magnetic 
field in order to bend the tracks of charged particles for momentum and charge sign 
determination. 

1.1.1 Tracking system 

The tracking system determines whether the particles are charged. In conjunction 
with a magnetic field, it measures the sign of the charge and the momentum of 
the particle. Photons may convert into an electron-positron pair and can in that 
case be detected in the tracking system. Moreover, charged kaon decays may be 
detected in a high-resolution tracking system through their characteristic "kink" 
topology: e.g. K='= — )• /U^f^ (64%) and K='= — tt^tt^ (21%). The charged parent 
(kaon) decays into a neutral daughter (not detected) and a charged daughter with 
the same sign. Therefore, the kaon identification process reduces to the finding of 
kinks in the tracking system. The kinematics of this kink topology allows to sepa- 
rate kaon decays from the main source of background kinks coming from charged 
pion decays 131 . 

1.1.2 Calorimeters 

Calorimeters detect neutral particles, measure the energy of particles, and deter- 
mine whether they have electromagnetic or hadronic interactions. EM and hadron 
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Figure 2: Comparison of the relative energy resolutions (as given by Eq. [T]) 
of the different EM calorimeters (left image) and hadronic calorimeters (right 
image) at the LHC experiments. The values of the parameters a, b and c were 
in all cases determined by fits to the data from beam tests and are given in the 
descriptions of the different experiments in Sects. |2.1[ |2.2| and |2.3[ In case 



of the ATLAS and CMS hadronic calorimeters the resolutions of the whole 
systems combining EM and hadronic calorimeters are shown. 



calorimetry at the LHC is described in detail in Refs. ||4l|5l. PID in a calorimeter 
is a destructive measurement. All particles except muons and neutrinos deposit all 
their energy in the calorimeter system by production of electromagnetic or hadronic 
showers. The relative resolution with which the deposited energy is measured is 
usually parametrised as: 




The first term takes into account the stochastic fluctuations and limits the low en- 
ergy performance. The second term is due to electronics noise. The third, constant 
term takes into account detector uniformities and errors in the calibration. This 
term limits the calorimeter performance at high energies. Fig. [2] shows a compari- 
son of the relative resolutions for the EM and hadronic calorimeters of the different 
LHC experiments. 

Calorimeters may be broadly classified into one of two types: sampling calorime- 
ters and homogeneous calorimeters. Sampling calorimeters consist of layers of 
a dense passive absorber material interleaved with active detector layers. In ho- 
mogeneous calorimeters on the other hand the absorber also acts as the detection 
medium. 

Photons, electrons and positrons deposit all their energy in the EM calorimeter. 
Their showers are indistinguishable, but an electrorj^can be identified by the exis- 
tence of a track in the tracking system that is associated to the shower. In this case 
the energy deposit must match the momentum measured in the tracking system. 
Hadrons on the other hand deposit most of their energy in the hadronic calorimeter 

'The term "electron" can sometimes refer to both "electron" and "positron" in this article. The 
usage should be clear from the context. The same is true for the muon and its anti-particle. 
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(part of it is also deposited in the EM calorimeter). However, the individual mem- 
bers of the families of charged and neutral hadrons can not be distinguished in a 
calorimeter. 

1.1.3 Muon system 

Muon systems at the LHC are described in detail in Ref. ||6l. The muon differs 
from the electron only by its mass, which is around a factor 200 larger. As a con- 
sequence, the critical energy Ec (the energy for which in a given material the rates 
of energy loss through ionization and bremsstrahlung are equal) is much larger for 
muons: it is around 400 GeV for muons on copper, while for electrons on cop- 
peij^it is only around 20MeV. As a consequence, muons do in general not produce 
electromagnetic showers and can thus be identified easily by their presence in the 
outermost detectors, as all other charged particles are absorbed in the calorimeter 
system. 

1.1.4 Other particles 

Neutrinos do in general not interact in a particle detector of the sort shown in Fig.[T| 
and therefore escape undetected. However, their presence can often be inferred by 
the momentum imbalance of the visible particles. In electron-positron colliders it 
is usually possible to reconstruct the neutrino momentum in all three dimensions 
and its energy. 

Quark flavor tagging identifies the flavor of the quark that is the origin of a 
jet. The most important example is B-tagging, the identification of beauty quarks. 
Hadrons containing beauty quarks can be identified because they have sufficient 
lifetime to travel a small distance before decaying. The observation of a secondary 
vertex a small distance away from the colUsion point indicates their presence. For 
this the information of a high-precision tracking system around the collision point 
is used. Such vertex tracking detectors are described in detail in Ref. [7J. Also 
Tau leptons, with a mean lifetime of 0.29 ps, fly a small distance (about 0.5 mm) 
before decaying. Again, this is typically seen as a secondary vertex, but without 
the observation of a jet. 

K^, A and A particles are collectively known as V*^ particles, due to their char- 
acteristic decay vertex, where an unobserved neutral strange particle decays into 
two observed charged daughter particles: e.g. — 7r+ vr^ and A — s-p vr^. V*^ 
particles can be identified from the kinematics of their positively and negatively 
charged decay products (see Figj3]l fE\. 

^For electrons an approximation for tiie critical energy is given by Ec — ^ MeV, where Z is 
the charge number of the material. 
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Figure 3: Armenteros-Podolanski plot fiom the ALICE experiment using data from 
^/s = 900 GeV proton collisions. The different Vq particles can be separated using 
the kinematics of their decay products. The orientation of the decay is described 
with respect to the momentum vector of the Vq. pf^ are the longitudinal momenta 
of the positively and negatively charged decay products with respect to the Vq 
particle's direction, qt represents the transverse component of the momentum of 
the positive decay product. 



1.2 PID by mass determination 

The three most important charged hadrons (pions, kaons and protons) and their an- 
tiparticles have identical interactions in an experimental setup as the one shown in 
Fig. [T] (charge deposit in the tracking system and hadronic shower in the calorime- 
ter). Moreover, they are all effectively stable. However, their identification can be 
crucial, in particular for the study of hadronic decays. The possible improvement 
in the signal-to-background ratio when using PID is demonstrated in Fig. [4| using 
the example of the (p — K+K^ decay. 

In B physics, the study of hadrons containing the beauty quark, different decay 
modes usually exist, and their individual properties can only be studied with an ef- 
ficient hadron identification, which improves the signal-to-background ratio (most 
tracks are pions from many sources). 

PID is equally important in heavy-ion physics. An example is the measure- 
ment of open charm (and open beauty), which allows to investigate the mecha- 
nisms for the production, propagation and hadronization of heavy quarks in the 
hot and dense medium produced in the collision of heavy ions. The most promis- 
ing channel is D'' — )• K~7r"'". It requires a very efficient PID, due to the small 
signal-to-background ratio. 

In order to identify any stable charged particle, including charged hadrons, it is 
necessary to determine its charge ze and its mass m. The charge sign is obtained 
from the curvature of the particle's track. Since the mass can not be measured 
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Figure 4: Demonstration of the power of PID by mass determination, using the 
example of the — K+K" decay measured with the LHCb RICH system (pre- 
Uminary data from t/s = 900 GeV p-p collisions |9]). The left image shows the 
invariant mass obtained from all combinations of pairs of tracks without PID. The 
right image shows how the (j) meson signal appears when tracks can be identified 
as kaons. 



directly, it has to be deduced from other variables. These are in general the mo- 
mentum p and the velocity fi = v/c, where one exploits the basic relationship 

P 

p = ■ymv — m = — — . (2) 

Here c is the speed of light in vacuum and 7 = (1 — /3^)~^/^ is the relativistic 
Lorentz factor. The resolution in the mass determination is 



my \P J V Py 
Since in most cases 7 » 1, the mass resolution is determined mainly by the accu- 
racy of the velocity measurement, rather than the momentum determination. 

The momentum is obtained by measuring the curvature of the track in the mag- 
netic field. The particle velocity is obtained by means of one of the following 
methods: 

1 . measurement of the energy deposit by ionization, 

2. time-of-flight (TOF) measurements, 

3. detection of Cherenkov radiation or 

4. detection of transition radiation. 



Each of these methods provide PID not only for charged hadrons, but also for 
charged leptons. The small obstacle of muons and pions not being well separated 
due to k, m.^ can luckily be circumnavigated, since muons can be easily iden- 
tified by other means. 
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Figure 5: Approximate minimum detector length required to achieve a K/vr separa- 
tion of > 3 a with three different PID techniques. For the energy loss technique 
we assume a gaseous detector. For the TOF technique, the detector length repre- 
sents the particle flight path over which the time-of-flight is measured. For the 
Cherenkov technique only the radiator thickness is given. The thicknesses of an 
expansion gap and of the readout chambers have to be added. 



The use of these methods is restricted to certain momentum ranges. For a given 
momentum range, the separation power can be used to quantify the usability of a 
technique. It defines the significance of the detector response R. If Ra and Rb 
are the mean values of such a quantity measured for particles of type A and B, 
respectively, and {(ta,b) is the average of the standard deviations of the measured 
distributions, then the separation power is given by 

n, = ^. (4) 

A summary of the momentum coverage and required detector lengths using the 
example of K/vr separation with the requirement > 3 is given in Fig. [5] 

Naturally, when choosing a PID technique, also other features have to be con- 
sidered besides the separation power. In practice, these often include luminosity 
and event rates, size and space requirements, accessibility, multiple scattering in 
the used materials, compatibility with other detector subsystems and geometrical 
coverage. 



2 Overview of the large LHC experiments 

The four main LHC experiments are the largest current collider experiments and 
integrate todays state-of-the-art detector technology, in particular with respect to 
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PID. This section gives an overview on the four experiments, with emphasis on the 
way the underlying physics program influences the experiment design with respect 
to PID. 

2.1 ATLAS and CMS 

ATLAS JTOl (A Toroidal Large AparatuS) and CMS lUl (Compact Muon Solenoid) 
are often called "general purpose" particle physics experiments, since they aim at 
uncovering any new phenomena appearing in proton-proton (p-p) collisions at the 
new energy domain that is now probed at the LHC. The main focus of the experi- 
ments is the investigation of the nature of the electroweak symmetry breaking and 
the search for the Higgs boson. The experimental study of the Higgs mechanism 
is expected to shed light on the mathematical consistency of the Standard Model 
(SM) of particle physics at energy scales above ^■l TeV. As a matter of fact, various 
alternatives to the SM foresee new symmetries, new forces, and new constituents. 
Furthermore, there are high hopes for discoveries that could lead towards a unified 
theory (super symmetry, extra dimensions, . . . ). 

2.1.1 Requirements 

The discovery and study of the Higgs boson is the benchmark process that influ- 
enced the design of the two experiments. Electrons, muons and photons are impor- 
tant components of its possible physics signatures. Assuming a low mass Higgs 
boson {rriH < 150GeV/c^), the predominant decay mode into hadrons is difficult 
to detect due to QCD backgrounds. In that case an important decay channel is 
H — )• 77. For masses above 130GeV, the most promising channel to study the 
properties of the Higgs boson is H ^ ZZ^*^ — )• 4/. It is called "gold-plated" in 
the particular case in which all four leptons are muons, due to the relative ease in 
detecting muons. The ATLAS and CMS detectors have been optimised to cover 
the whole spectrum of possible Higgs particle signatures. In summary, these are 
the requirements: 

1 . large acceptance in pseudo-rapiditjj^ and almost full coverage in azimuthal 
angle (the angle around the beam direction); 

2. good identification capabilities for isolated high transverse momentunj^pho- 
tons and electrons; 

3. good muon ID and momentum resolution over a wide range of momenta 
and angles. At highest momenta (1 TeV) a transverse momentum resolution 
(Tpj./pT of the order 10% is required. 

^The pseudo-rapidity is defined as = — ln[ian(0/2)], where Q is tiie polar angle between the 
charged particle direction and the beam axis. 

''The transverse component of the momentum is the one in the plane that is perpendicular to the 
beam direction. 
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To maximize the integrated luminosity for the rare processes that are the main 
interest of the experiments, proton bunches will be brought to collision every 25 ns. 
Strong focusing of the beams helps to increase the luminosity to lO^^cm^^s^^. 
The consequences are a large number of overlapping events per proton bunch cross- 
ing and extreme particle rates. For the experiments these pose major challenges. 

2.1.2 Setup 

The setup of ATLAS and CMS is in general quite similar, following the "tradi- 
tional" setup as described in Sect. |1.1[ However, in the implementation of some 
of the components quite different choices were made. The main similarities and 
some differences can be seen in Tab. [T] Schematic images of the two detectors are 
shown in Figs. [6] and [7] 

2.1.3 Tracking and muon systems 

As the innermost component, the tracking systems of ATLAS and CMS are em- 
bedded in solenoidal magnetic fields of 2 T and 4 T, respectively. In both cases the 
tracking systems consist of silicon pixel and strip detectors. ATLAS includes also 
a Transition Radiation Tracker (TRT) based on straw tubes, which provides also 
electron ID [2|. 

The global detector dimensions are defined by the large muon spectrometers, 
which are designed to measure muon momenta with extremely high accuracy ||6l. 
While the muon system of ATLAS is designed to work independently of the inner 
detector, in CMS the information from the tracking system is in general combined 
with that from the muon system. 

The ATLAS muon system uses eight instrumented air-core toroid coils, provid- 
ing a magnetic field mostly orthogonal to the muon trajectories, while minimising 
the degradation of resolution due to multiple scattering (T0\. For high-precision 
tracking, the magnets are instrumented with Monitored Drift Tubes (MDT) and, 
at large pseudo-rapidities. Cathode Strip Chambers (CSC). By measuring muon 
tracks with a resolution <50 /im, a standalone transverse momentum resolution 
(t{pt)/pt ~ 3% is achieved for = lOOGeV/c, while at px = ITeV/c it is 
around 10%. Only below 200GeV/c the combination of the information from the 
muon and tracking systems helps, keeping the resolution below 4%. As a separate 
trigger system and for second coordinate measurement. Resistive Plate Chambers 
(RPCs) are installed in the barrel, while Thin Gap Chambers (TGCs) are placed in 
the end-caps, where particle rates are higher. 

The CMS design relies on the high bending power (12 Tm) and momentum 
resolution of the tracking system, and uses an iron yoke to increase its magnetic 
field nil. With the field parallel to the LHC beam axis, the muon tracks are bent 
in the transverse plane. The iron yoke is instrumented with aluminum Drift Tubes 
(DT) in the barrel and CSCs in the end-cap region. Due to the iron yoke, the 
momentum resolution of the CMS muon system is dominated by the multiple scat- 
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Figure 6: Perspective view of the ATLAS detector flO]. The dimensions are 25 m 
in height and 44 m in length, the overall weight of the detector is approximately 
7000 tonnes. 



Table 1 : Overview on the technologies chosen for the ATLAS and CMS tracking 
systems, EM and hadronic calorimeters and muon systems. Their acceptances in 
pseudo-rapidity r] are given as well. Abbreviations are explained in the text. 



Detector 


technology 
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component 
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Si pixel detector (3 layers) 








Si strip detector 
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(Tracking) 


MDT 
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Muon 
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|r/| < 1.05 


hi < 1-6 


(Trigger) 


TGC 1 


1.05 < I77I < 2.7 
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Figure 7: Perspective view of the CMS detector ifTTI . The dimensions are 14.6 m 
in height and 21.6 m in length. The overall weight is approximately 12 500 tonnes. 

tering. As a consequence, the requirements on spatial resolution are somewhat 
looser (~70 ^m) than in ATLAS. The standalone muon momentum resolution is 
(^{pt)/pt = 9% iorpT < 200GeV/c and 15 to 40% aipT = ITeV/c, depending 
on \r]\. Including the tracking system improves the result by an order of magnitude 
for low momenta. At 1 TeV the contribution of both measurements together yield a 
momentum resolution of about 5%. The DT and CSC subsystems can each trigger 
on muons with large transverse momentum. However, for the full LHC luminos- 
ity, faster trigger chambers are needed to associate the detected muons to the right 
crossing of proton bunches in the LHC. RPCs are used throughout the whole CMS 
muon system for that purpose. 

2.1.4 Particle identification 

Electrons, hadrons and neutral particles are identified in the calorimeter systems, 
while muons are identified in the large muon systems. The designs of the EM 
calorimeters and muon systems have been guided by the benchmark processes of 
Higgs boson decays into two photons or into leptons. 

The ATLAS EM calorimeter is of sampling type and consists of liquid Ar- 
gon (LAr) as detection medium and lead (Pb) absorber plates operated in a cryo- 
stat at 87 K. The calorimeter depth varies from 22 to 38 Xq, depending on the 
pseudo-rapidity range. In beam tests [ 10] the relative energy resolution after sub- 
traction of the noise term was found to be {(TE/Ef = (0.1/y£^^GeV))^+0.0072. 
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The CMS EM calorimeter on the other hand is a homogeneous calorimeter made 
from lead tungstate (PbW04) crystals with a length that corresponds to around 
25 Xq. The relative energy resolution was measured in beam testj^as {ge/EY = 
(0.028/v^-E; (GeV))2 + (0.125/S( GeV))^ + O.OOS^. 

The ATLAS hadronic calorimeter is made from a Barrel and two end-cap 
modules. The barrel part is made from steel and plastic scintillator tiles. The 
end-caps are made from copper plates and use LAr as active medium. In beam 
tests of the combined EM and hadronic calorimeter systems the relative energy 



resolution for pions was found to be {aE/Ef = {Q.52/yjE (GeV))^ + 
(0.016/£;(GeV))2 + 0.03^ in the barreQ and {as/Ef « (0.84/y/^ (GeV))^ in 
the end-caps. The CMS hadronic calorimeter is made from brass and plastic scintil- 
lator tiles. The resulting energy resolution ||T2|| for the combined system of EM and 
hadronic calorimeters and for pions is {aE/Ef = {1.12/ ^JE (GeV))^ + 0.036^. 

The silicon detectors of the ATLAS and CMS tracking systems offer the pos- 
sibility of hadron ID at low momenta (a few hundred MeV/c) via ionization mea- 
surements. For electron ID at momenta up to 25 GeV/c ATLAS also takes into 
account the information from the TRT, namely large energy deposits by electrons 
due to transition radiation (for details see Ref. 111). 



2.2 ALICE 

ALICE [ 13] (A Large Ion Collider Experiment) is the dedicated heavy-ion experi- 
ment at the LHC. The LHC can collide lead (Pb) nuclei at center-of mass energies 
of sJsnn = 2.76 and 5.5 TeV. This leap to up to 28 times beyond what is presently 
accessible will open up a new regime in the experimental study of nuclear mat- 
ter. The aim of ALICE is to study the physics of strongly interacting matter at the 
resulting extreme energy densities and to study the possible formation of a quark- 
gluon plasma. Many observables meant to shed light on the evolution of the quark- 
gluon plasma depend on FID. The most natural example is particle spectra, from 
which the freeze-out parameters (such as the kinetic and chemical freeze-out tem- 
perature and the collective flow velocity) can be extracted. Besides pions, kaons 
and protons, also resonances such as the (j) meson can be used for such analysis. 
Actually, any modification in the mass and width of the meson could be a sig- 
nature for partial chiral symmetry restoration, and thus, creation of a quark-gluon 
plasma. The resonance is typically identified via its hadronic decay channel 
(f) — )• K~ , which makes a good hadron ID important. Leptonic channels are 
very important as well and thus extremely good lepton ID is needed. 

ALICE also studies p-p collisions at y/s = 7 and 14 TeV. These data are im- 
portant as a reference for the heavy-ion data. Moreover, due to the very low mo- 

^The measurement was done for electrons, requiring their impact point to lie in the center of the 
3x3 crystals used to evaluate the energy deposit. Without this requirement the relative resolution is 
slightly worse but still meeting the design goal of better than 0.5% for E > 100 GeV [ 1 1 1. 

*For the barrel hadron calorimeter standalone the relative energy resolution for pions is 
{(JE/Ef = (0.564/7S (GeV))2 + 0.055^ 
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mentum cut-off (~0. IGeV/c), ALICE can contribute in physics areas where it 
complements the other LHC experiments. To name an example, the measurement 
of charm and beauty cross sections is possible down to very low momenta. More- 
over, proton physics at high multiplicities is easily accessible via the multiplicity 
trigger from the Silicon Pixel Detector (SPD). 

2.2.1 Requirements 

As compared to p-p collisions, the charged track multiplicities in Pb-Pb collisions 
are extraordinary. At the time of the ALICE Technical Proposal, charged particle 
densities of up to dNch/drj = 8000 were considered for the LHC center-of-mass 
energy |[T4l [T5l [T6l . Including secondaries, this would amount to 20 000 tracks in 
one interaction in the relevant acceptance region. However, first measurements ifTTl 
at ^/~sn~n = 2.76 TeV yield dNch/dr] ^ 1600. When compared to p-p collisions 
at the LHC, the produced particles have relatively low momenta: i.e. 99% have a 
momentum below about 1 GeV/c. Hadrons, electrons, muons and photons may be 
used as probes in order to explore the strongly interacting matter that is produced 
in the Pb-Pb collisions. Those probes are inspected by dedicated PID detectors. 
In some cases these detectors may have very limited geometrical acceptance; for 
some actually not even a coverage of 2 vr in azimuthal angle is necessary. This 
is an important difference when compared to general purpose experiments. Ne- 
vertheless, large geometrical acceptance remains important for the main tracking 
devices. In summary, these are the requirements for the ALICE experiment: 

1 . reliable operation in an environment of very large charged track multiplici- 
ties; 

2. precision tracking capabilities at very low momenta (lOOMeV/c), but also 
up to lOOGeV/c; 

3. low material budget; 

4. low magnetic field (0.2 < B < 0.5 T) in order to be able to track low 
momentum particles; 

5. good hadron ID for momenta up to a few GeV/c and electron ID up to 
lOGeV/c in the central barrel; 

6. good muon ID (in the forward region). 

At LHC energies the total cross-section for Pb-Pb collisions will be of the order 
of 8 b. At the maximum expected luminosity of 10^^cm~^s^^ this corresponds 
to an interaction rate of 8 kHz. Only a fraction of these events will be central 
collisions (collisions with small impact parameter) and the aim is to be able to 
trigger on at least 200 of such central events per second. 
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Figure 8: Perspective view of the ALICE detector fTT]. The dimensions are 16m 
in height and 26 m in length. The overall weight is approximately 10 000 tonnes. 



For p-p collisions the goal is to be able to read out the experiment at rates of 
at least 1.4 kHz. Since the specific features of the used detectors (TPC drift time 
of almost 100 /is, see Sect. |3.7| ) makes it a rather slow detector (at least compared 
to the three other large LHC experiments), the luminosity for p-p collisions has 
to be limited [15|: at ~5 x 10^" cm^^s^^, corresponding to an interaction rate of 
~200kHz, the integrated luminosity can be maximised for rare processes. With 
lower luminosity (10^^ cm~^s~^) statistics for large cross section observables can 
be collected and global event properties can be studied at optimum detector perfor- 
mance. 



2.2.2 Setup 

The ALICE experimental setup is shown in Fig. [8] ALICE consists of a central 
barrel part for the measurement of hadrons, leptons and photons, and a forward 
muon spectrometer. The central part is embedded in a large solenoid magnet that 
is reused from the L3 experiment at the LEP collider at CERN. 

A Time Projection Chamber (TPC) was chosen as the main tracking device in 
the central barrel, as this is a very high granularity detector capable of satisfying the 
requirements given in the previous section. The TPC surrounds the Inner Tracking 
System (ITS) which is optimised for the determination of primary and secondary 
vertices and for precision tracking of low-momentum particles. As the innermost 
layer, the Silicon Pixel Detector (SPD) has a key role in the determination of the 
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vertex position. A unique feature of the SPD is its capability to generate a prompt 
trigger based on a programmable, fast online algorithm. As second and third lay- 
ers, a Silicon Drift Detector (SDD) and a Silicon Strip Detector (SSD) complete the 
ITS, with two detection planes each. Finally, on the outside of the TPC, the Transi- 
tion Radiation Detector (TRD) contributes to the tracking capability and provides 
electron ID El. 



2.2.3 Particle identification 

The strategy for PID in the barrel of ALICE is described in detail in Ref. fTE\. 
In a first step, all charged particles from the collision are tracked in the ITS, TPC 
and TRD using parallel Kalman filtering. The PID procedure is then applied to all 
reconstructed tracks that have been successfully associated to a signal in one of the 
PID detectors. 

The TPC, SSD and SDD each provide PID via ionization measurements. The 



TPC is described in more detail in Sect. 3.7 The TRD is designed for electron 



ID and is described in Ref. ||2l. The time-of-flight (TOF) array provides charged 
hadron ID and is described in Sect. 14.51 

Surrounding the TOF detector, there are three single-arm detectors inside the 
ALICE solenoid magnet: the Photon Spectrometer (PHOS), the Electro-Magnetic 
CALorimeter (EMCAL) and an array of RICH counters optimised for High-Momentum 



Particle IDentification (HMPID, see Sect. 5.3 1. PHOS is a homogeneous EM calo 



rimeter based on lead tungstate crystals, similar to the ones used by CMS, read 
out using Avalanche Photodiodes (APD). It covers 100 deg in azimuthal angle and 
|?7| < 0.12 in pseudo-rapidity. Its thickness corresponds to 20 Xq. The crystals 
are kept at a temperature of 248 K to minimize the contribution of noise to the 
energy resolution, optimizing the response for low energies. In a beam test the 



relative energy resolution was found to be {cte/EY = (0.033/i/i? (GeV))^ + 
(0.018/^; (GeV))2 + 0.011^. The EMCAL was proposed as a late addition to the 
ALICE setup, with the primary goal to improve the capabilities of ALICE for jet 
measurements. The main design criterion was to provide as much EM calorimetry 
coverage as possible within the constraints of the existing ALICE detector sys- 
tems. The EMCAL is a sampling calorimeter made from lead absorber plates and 
scintillators and is positioned opposite in azimuth to the PHOS, covering 107 deg 
in azimuthal angle and \ri\ < 0.7 in pseudo-rapidity. The relative energy reso- 
lution of the EMCAL was measured in a test beam and can be parameterised as 



{(TE/Ef = (0.113/^^ (GeV))2 + 0.01682. 

The goal of the ALICE muon spectrometer is to study vector mesons contain- 
ing heavy quarks, such as J/^', and members of the T family via the muonic 
channel. The muon spectrometer covers only the limited pseudo-rapidity interval 
2.5 < ?7 < 4. Closest to the interaction region there is a front absorber to remove 
hadrons and photons emerging from the collision. Five pairs of high-granularity 
detector planes form the tracking system within the field of a large dipole magnet. 
Beyond the magnet is a muon filter (a 120 cm thick iron wall), followed by two 
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pairs of trigger chamber planes (RPCs). 
2.3 LHCb 

LHCb lITSl is dedicated to heavy flavor physics. One particular aim is to look 
at evidence of new physics in CP violation and rare decays of beauty and charm 
hadrons. The level of CP violation in the SM can not explain the absence of anti- 
matter in our universe. A new source of CP violation is needed to understand this 
matter-antimatter asymmetry, implying new physics. Particles associated with new 
physics could manifest themselves indirectly in beauty or charm meson decays and 
produce contributions that change the expectations of CP violation phases. They 
may also generate decay modes forbidden in the SM. 

2.3.1 Requirements 

A large bb production cross section of the order ~500 ^b is expected for p-p col- 
lisions at ^/s = 14TeV. At high energies, the bb pairs are predominantly pro- 
duced in a forward or backward cone. Separating pions from kaons in selected 
B hadron decays is fundamental to the physics goals of LHCb. A good example 
is the channel Bg — D^K^, which has to be separated from the background 
from Bg — )• DgiT^, which is about 15 times more abundant. Moreover, other fi- 
nal states containing electrons, muons and neutral particles (photons, neutral pions 
and ry's) have to be distinguishable. The requirements for the LHCb detector can 
be summarised in this way: 

1. geometrical acceptance in one forward region (1.9 < t] < 4.9); 

2. good hadron ID; 

3. good momentum and vertex resolutions and 

4. an efficient and flexible trigger system. 

In order to maximize the probability of a single interaction per beam crossing, 
the luminosity in the LHCb interaction region may be limited to 2 to 5 x 10'^^ cm^s^. 
In these conditions, one year of LHCb running (~10^ s) corresponds to 2 fb^^ of 
integrated luminosity and about 10^^ bb pairs produced in the region covered by 
the spectrometer. 

2.3.2 Setup 

Unlike ATLAS and CMS, LHCb does not have a cylindrical geometry. Rather, it 
is laid out horizontally along the beam line, as shown in Fig. [9] The tracking sys- 
tem consists of the VErtex LOcator (VELO) and four planar tracking stations: the 
Tracker Turicensis (TT) upstream of the 4Tm dipole magnet, and stations T1-T3 
downstream of the magnet. VELO and TT use silicon strip detectors. In T1-T3, 
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Figure 9: Schematic view of the LHCb detector ifTSl . 



silicon strips are used in the region close to the beam pipe, whereas strawtubes 
are employed in the outer regions. The VELO makes possible a reconstruction of 
primary vertices with 10 /xm (60 ^m) precision in the transverse (longitudinal) di- 
rection. In this way the displaced secondary vertices, which are a distinctive feature 
of beauty and charm hadron decays, may be identified. The overall performance 
of the tracking system enables the reconstruction of the invariant mass of beauty 
mesons with resolution a-m ~ 15 to 20 MeV/c^, depending on the channel. 

2.3.3 Particle identification 

LHCb in general looks like a slice out of a "traditional" experiment as described in 
Sect. |1.1[ apart from the two RICH detectors providing hadron ID. The RICH de- 
tectors are described in more detail in Sect. 15.41 An EM calorimeter and a hadron 
calorimeter provide the identification of electrons, hadrons and neutral particles 
(photons and ttq) as well as the measurement of their energies and positions. The 
EM calorimeter is a rectangular wall constructed out of lead plates and scintilla- 
tor tiles. The total thickness corresponds to 25 Xq. In a beam test it was found 
that the relative energy resolution follows (aE/Ef = (0.094/i/£' (GeV))^ + 
(0.145/-E (GeV))^ + 0.0083^. The hadronic calorimeter consists of iron and scin- 
tillator tiles with a relative energy resolution of {(Te/E)"^ = (0.69/y^£' (GeV))^ + 
0.09^, measured with a prototype in a beam test. Finally, the muon system is de- 
signed to provide a fast trigger on high momentum muons as well as offline muon 
indentification for the reconstruction of muonic final states and beauty flavor tag- 
ging. It consists of five stations (M1-M5) equipped mainly with Multi Wire Pro- 
portional Chambers (MWPCs). For the innermost region of station Ml, which has 
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the highest occupancy, Gas Electron MultipUers (GEMs) |[T9l are used. For a muon 
PID efficiency of 90% the misidentification rate is only about 1.5%. 



3 Ionization measurements 

Ionization of matter by charged particles is the primary mechanism underlying 
most detector technologies. The characteristics of this process, along with the 
momentum measurement, can be used to identify particles. 

3.1 Energy loss and ionization 

When a fast charged particle passes through matter, it undergoes a series of inelastic 
Coulomb collisions with the atomic electrons of the material. As a result, the 
atoms end up in excited or ionised states, while the particle loses small fractions 
of its kinetic energy. The average energy loss per unit path length {AE/Ax) is 
transformed into the average number of electron/ion pairs (or electron/hole pairs 
for semiconductors) {Nj) that are produced along the length x along the particle's 
trajectory ll20ll : 



where W is the average energy spent for the creation of one electron/ion (elec- 
tron/hole) pair. W exceeds the ionization energy Ei (or the band gap energy Eg 
for a solid) of the material, because some fraction of the energy loss is dissipated 
by excitation, which does not produce free charge carriers. Typical values of W 
lie around 30 eV for gases, being constant for incident particles with relativistic 
velocities (/3 1), but increasing for low velocities. For semiconductors the W 
values are roughly proportional to the band gap energy: 



and are much lower than for gases 1211 : e.g. on average 3.6 eV in silicon and 
2.85 eV in germanium. Consequently, the ionization yield in semiconductor detec- 
tors is much larger than in gaseous devices. 

The interactions of the charged particle with the atomic electrons can be mod- 
eled in terms of two components: primary and secondary interactions. In primary 
interactions direct processes between the charged particle and atomic electrons 
lead to excitation or ionization of atoms, while secondary processes involve subse- 
quent interactions. The primary interactions can be characterised by the Rutherford 
cross-section (with the energy dependence Aa /dE oc £'~^) for energies above the 
highest atomic binding energy, where the atomic structure can be ignored. In this 
case the particle undergoes elastic scattering on the atomic electrons as if they 
were free. According to the steeply falling Rutherford spectrum most of the pri- 
mary electrons emitted in such colhsions have low energy. However, a significant 




(5) 



W = 2.^Eg + {).QdW 



(6) 
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probability for producing primary electrons with energies up to the kinematic limit 
for the energy transfer Emax exists. Emax is given by 



where rrif, is the electron mass, x = nie/m and m is the mass of the incident 
particle. In such collisions, characterised by a very small impact parameter, the 
energy transfered to the electron will be larger than Ej (or Eg) and the resulting 6- 
rays or knock-on electrons produce additional ionization in secondary interactions. 
5-rays can even leave the sensitive volume of the detector, but a magnetic field may 
force them to curl up close to the primary charged particle's track. In this case they 
will contribute to a measurement of the deposited energy. 

In colhsions with large impact parameter the atomic electrons receive much 
less energy, which is used for excitation without the creation of free charges. How- 
ever, in gases tertiary ionization by collisions of an atom in an excited state with 
other atoms may be important (Penning ionization) [.22il . 

3.2 Velocity dependence 

The first calculation for the average energy loss per unit track length based on the 
quantum mechanical principles of the scattering theory was introduced by Hans 
Bethe in 1930 and 1932 [24*, "^S"]. The well known Bethe-Bloch formula is mod- 
ified to yield the restricted (average) energy loss by neglecting higher energy 5- 
electrons through the introduction of an upper limit for the energy transfer in a 
single collision Ecut li20l : 



Here ze is the charge of the incident particle and I is the effective excitation en- 
ergy of the absorber materia^ measured in eV. 5 is the density effect correction to 
the ionization energy loss, which was calculated for the first time by E. Fermi in 
1939 L26|. It is much larger for liquids or solids than for gases and is usually com- 
puted using a parameterization by R. M. Sternheimer [23 1. Typical values for E^ut 
depend on the strength of the magnetic field and lie in the range 10 to lOOkeV. 

Eq. [8] is valid for electrons and also for heavier particles. In the low energy 
region (i.e. < 0.5) the average energy loss decreases like l//3^ and reaches 
a broad minimum of ionisation (MIP region) at f]j ^ 4. For larger values of 
the average energy loss begins to rise roughly proportional to log(/37) (relativistic 
rise) with a strength defined by the mean excitation energy I. The rise is reduced 
at higher momenta by the density effect correction 6. A remaining relativistic rise 

'For elements the excitation energy (in eV) can be approximately calculated as J ~ 13.5 Z for 
Z < 14 and I ^ 10 Z for Z > 14, where Z is the atomic number of the absorber material. A 
graphical presentation of measured values is given in Ref. 1231 . 
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would be due to (rare) large energy transfers to a few electrons. Since these events 
are removed in Eq. [8] by the introduction of Ecut, the restricted average energy 
loss approaches a constant value (the Fermi plateau) as /Sj increases. For solids, 
due to a stronger density effect, the Fermi plateau is only a few percent above the 
minimum. 

3.3 Straggling functions 

For a given particle, the actual value of the energy loss A over a given path length 
X is governed by statistical fluctuations which occur in the number of collisions 
(a Poissonian distribution) and in the energy transferred in each collision (a dis- 
tribution ~1/E^). The resulting distribution of the energy loss F{x, A) is called 
straggling function. One particular way to calculate this distribution using specific 
assumptions was introduced by L. Landau in 1944 [27]. However, in particle phy- 
sics often the name Landau junction is used generically to refer to all straggling 
functions. 

The ionization distribution (the distribution of the number of electron/ion or 
electron/hole pairs Nj along the particle path length x) G{x, Nj) can be obtained 
from F{x, A) by using the relation Nj = A/W, assuming that the energy loss is 
completely deposited in the material (the sensitive volume of the detector). Also 
G{x,Nj) is called a straggling function. 

The discussion of straggling functions is generally divided into two cases: thick 
absorbers and thin absorbers. In a thick absorber, with a thickness sufficient to ab- 
sorb an amount of energy comparable to the particles initial kinetic energy, the 
straggling function can be approximated by a Gaussian distribution. For thin ab- 
sorbers, where only a fraction of the kinetic energy is lost through ionizing colli- 
sions, the straggling function is always influenced by the possibility of large energy 
transfers in single collisions. These add a long tail (the Landau tail) to the high en- 
ergy side, resulting in an asymmetric form with a mean value significantly higher 
than the most probable value. 

3.4 PID using ionization measurements 

Measurements of the deposited energy A can be used for PID. Gaseous or sohd 
state counter^ provide signals with pulse height R proportional to the number of 
electrons liberated in the ionization process along the track length x inside the 
sensitive volume, and thus proportional to A. To limit deterioration of the resolu- 
tion due to energy loss fluctuations, in general pulse height measurements are 
performed along the particle track, either in many consecutive thin detectors or in a 
large number of samples along the particle track in the same detector volume. For 
sufficiently large Nji the shape of the obtained pulse height distribution approaches 
the one of G{x,Nj). The mean value of the distribution is not a good estimator 

^In a typical particle physics experiment the counters used for ionization measurements are the 
same devices that also provide the spatial coordinates for the momentum measurement. 
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Figure 10: Typical curves of the ionization signal as a function of particle mo- 
mentum for a number of known charged particles. A parameterization like the one 
suggested in Ref. ||20|| was used to calculate the curves. 



for the energy deposit A, since it is quite sensitive to the number of counts in the 
tail, which limits the resolution. A better estimator is derived from the peak of the 
distribution, which is usually approximated by the truncated mean {R)a, which is 
defined as the average over the M lowest values of the pulse height measurements 

Ri- 

1 

i=l 

Here Ri < -Rj+i for z = 1, . . . , n — 1 and M is an integer M = uNr. a is a 
fraction typically between 0.5 and 0.85. Values of {R)a follow an almost perfect 
Gaussian distribution with variance aE- 

From a theoretical point of view it should be possible to achieve better results 
by using all the information available, i.e. the shape of G{x,Nj). A maximum 
likelihood fit to all Nji ionization measurements along a track or tabulated energy 
loss distributions may be used. However, in practice they show a performance that 
is not significantly better ll20l l28ll . Thus usually the truncated mean method is 
chosen because of its simplicity. 

In practice, only relative values of the deposited energy A = A(m, (3) are 
needed to distinguish between different particle species. Because Eq. [8] is not a 
monotonic function of /3, it is not possible to combine it with Eq. |2] in order to 
obtain a form m = m{A,p), where p is the momentum of the particle. However, 
the PID capability becomes obvious by plotting ionization curves A(m, (3) for a 
number of charged particles. A description of measured ionization curves based on 
five parameters is given in Ref. [20 1 and can be used for this purpose. An example 



is shown in Fig. 10 By simultaneously measuring p and A for any particle of un- 
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Figure 11: Typical separation power achievable with ionization measurements in a 
gaseous detector. The ionization curves from Fig. 10 were used together with an 
assumed energy resolution of 5%. 



known mass, a point can be drawn on this diagram. The particle is identified when 
this point can be associated with only one of the curves within the measurement 
errors. Even in areas where the bands are close, statistical PID methods may be 
applied. 



3.5 Energy resolution and separation power 

The resolution in the ionization measurement {energy resolution) (Je is given by 
the variance of the Gaussian distribution of the truncated mean values and is in 
general proportional to the energy deposit: cr^; oc A. Using Eq. |4j the separation 
for two particles A and B with different masses but the same momentum can be 
calculated as 

^''E = ~u ^ ■ ^^^^ 

\^a,b) 

The average of the two resolutions {(ta,b) = {(^E^A + o'£;,b)/2 is used. For the 
purpose of qualitatively discussing the PID capabilities this approximation is suffi- 
cient. For alternative ways to describe the resolution of the PID, refer to Ref. [28,1 . 



In Fig. 1 1 typical particle separations are shown for a gaseous detector with 
an energy resolution of 5%. As expected, hadron identification works well in the 
low momentum region. In the minimum ionization region, where hadrons carry 
momenta around a few GeV/c and the ionization curves are very close, it is likely 
that the method fails to discriminate the particles. In the region of the relativistic 
rise moderate identification is possible on a statistical basis. 
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3.6 Errors affecting the resolution 

The operational regime of the particle detector used for the ionization measure- 
ments must be chosen such that the measured pulse heights R are exactly propor- 
tional to the number of electron/ion (electron/hole) pairs Nj created by the energy 
deposit A. The same is true for the signal processing chain. Nevertheless, some 
effects change the apparent energy deposit and thus limit the energy resolution of 
the device. They have to be considered on top of the fundamental limit that is given 
by the statistics of the primary ionization. 

• The ionization signal amplitude is determined by the detector (e.g. gas am- 
plification) and by the electronics. An optimal performance can be achieved 
only if an energy calibration is carried out in order to determine the absolute 
gain of each channel with high precision. 

• Overlapping tracks must be eliminated with some safety margin. This can 
introduce some limit for a high track density environment like heavy-ion col- 
lisions or jets, because the number of ionisation measurements Nr entering 
the truncated mean calculation (Eq. [9]) can be considerably reduced by this 
effect. 

• Track independent effects have to be kept under control. For a gaseous de- 
vice these include the gas density (pressure and temperature) and gas mix- 
ture. Electronics effects like noise and baseline shifts have to be taken care 
of as well. 

• On the level of each track, more effects influence the apparent ionization: the 
detector geometry together with the track orientation influence the amount of 
ionization charge per sampling layer and can be corrected for. In a gaseous 
device, attachment of the drifting electrons decreases the ionization signal 
and requires a correction, which depends on the drift length z through a 
constant b: R{z) = Re~^^. A similar decrease of the signal as a function of 
the drift length can be due to diffusion of the drifting primary electron cloud, 
and can be corrected for as well. 

For gaseous detectors, these effects are described in detail in Ref. |[20l . An 
empirical formula exists to estimate the dependence of the energy resolution on 
the number of measurements Nr, on the thickness of the sampling layers x, and 
on the gas pressure P 131] |33 : 

cr£; = 0.41 7V-°-^3(3;P)-°-32 _ (11) 

The formula was obtained by fits to measured data and includes optimizations such 
as the truncated mean method and the ones just described. If the ionization mea- 
surements along a track were independent, and if no other error sources like elec- 
tronics noise existed, the resolution would scale as Nn^'^. Moreover, the formula 
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shows that for a fixed total length of a detector x Np;, one obtains a better reso- 
lution for a finer sampling, provided that the ionization is sufficient in each sam- 
pling layer. A comparison of the expected resolutions calculated with Eq. 1 1 and 
measured values is given in Ref. [33]. Summaries and examples of performances 
achieved with gaseous tracking detectors like Drift, Jet or Time Projection Cham- 
bers are given in Refs. ll20l l32l . Typical values are 4.5 to 7.5%. By increasing 
the gas pressure, the resolution can be significantly improved, e.g. for the PEP 
TPC an energy resolution of below 3% was achieved at 8.5 bar. However, at higher 
pressure the relativistic rise is reduced due to the density effect. 



3.7 ALICE Time Projection Chamber 

The ALICE TPC (Fig. [12]) is the largest gaseous Time Projection Chamber 

built so far. A TPC is a remarkable example of a detector performing precise track- 
ing and PID through measurements of the energy deposited through ionization. 
Since it is basically just a large volume filled with a gas, a TPC offers a maximum 
of active volume with a minimum of radiation length. In a collider experiment 
like ALICE, the field cage of a TPC is typically divided into two halves separated 
by a planar central electrode made by a thin membrane. The electrons produced 
by charged particles crossing the field cage drift towards the two end-caps, where 
readout detectors are mounted. In general, a gated wire grid is installed directly in 
front of the readout chambers. When the system is triggered, the gate opens, allow- 
ing the passage of the ionization electrons, which create charge avalanches in the 
readout detectors. Closing the gate assures that the ions created in the avalanche 
process do not enter the drift volume. Measuring the signals induced on adjacent 
readout pads makes possible an accurate determination of the position by evalu- 
ating the center of gravity. Together with accurate measurements of the arrival 
time (relative to some external reference such as the collision time of the beams) 
the complete trajectory of all charged particles traversing the TPC can be deter- 
mined. Due to the little amount of dead volume and the easy pattern recognition 
(continuous tracks), TPCs are the best tracking devices for high multiplicity envi- 
ronments, e.g. heavy-ion experiments. 



3.7.1 Field cage and readout chambers 

The ALICE TPC is a large-volume TPC with overall "conventional" lay-out but 
with nearly all other design parameters pushed towards the limits. A comprehen- 
sive description can be found in Ref. [30|. The TPC field cage is a hollow cylinder 
whose axis is aligned with the LHC beam axis and is parallel to the magnetic field. 
The active volume has an inner radius of about 85 cm, an outer radius of about 
250 cm, and an overall length along the beam direction of 500 cm. The central 
electrode is charged to 100 kV and provides, together with a voltage dividing net- 
work at the surface of the outer and inner cylinder, a precise axial electric field of 
400 V/cm. To ensure low diffusion of the drifting electrons and a large ion mobility. 
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Figure 12: 3D view of the TPC field cage ll30ll . The high voltage electrode is 
located at the center of the 5 m long drift volume. The two endplates are divided 
into 18 sectors holding two readout chambers each. 

Ne was chosen as the main component of the counting gas. The gas mixture with 
10 % of CO2 as quencher makes mandatory a very good temperature homogenisa- 
tion in the drift volume (AT < 0.1 K), in order to guarantee a homogeneous drift 
velocity. It is operated at atmospheric pressure. 

Wire chambers with cathode pads are used for the readout of the TPC. An 
optimization of the design for the expected high track density environment impUed 
rather small pad sizes (4 x 7.5 mm^ in the innermost region). As a consequence, the 
number of pads is large (~560 000; for comparison: the ALEPH TPC had 41 000 
pads), and the readout chambers have to be operated at rather high gas gains near 

3.7.2 Front-end electronics 

Also for the design of the front-end electronics the key input was the operation 
in a high track density environment [30|. A signal occupancy as high as 50% 
can be expected for central Pb-Pb collisions. Due to the large raw data volume 
(750MB/event) zero suppression is implemented in the front-end electronics, thus 
reducing the event sizes to about 300 kB for p-p collisions (no overlying events) 
and to about 12 MB for Pb-Pb collisions. Signal tail cancellation and baseline 
correction are performed before the zero suppression in order to preserve the full 
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Figure 13: Measured ionization signals of charged particles as a function of the 
track rigidity (particle momentum divided by charge number). 11x10® events from 
a data sample recorded with y/s = 7 TeV p-p collisions provided by the LHC were 
analysed. The lines correspond to a parameterization of the Bethe-Bloch curve as 
described in Ref. |20|. Also heavier nuclei like deuterons and tritium and their 
anti-particles are found. 



resolution on the signal features (charge and arrival time). A large dynamic range 
ensures that the ionization signal of particles can be determined with precision 
from very low up to high momenta. 



3.7.3 PID performance and outlook 

A very good energy calibration is fundamental for an optimal PID performance (see 



Sect. 3.6 1. In the ALICE TPC, the amplitude of the ionization signal is determined 



by the gas amplification and the electronics. The absolute gain of each channel is 
obtained with high precision by periodically releasing radioactive krypton (fgKr) 
into the TPC gas system and measuring the well known decay pattern (for more 
details see ll30ll ). 

Inside the TPC, the ionization strength of all tracks is sampled on up to 159 
pad rows. The energy resolution is 5% for tracks with the maximum number of 



signals. Fig. 13 shows the ionization signals from 11 x 10® events from a data 
sample recorded with y/s = 7 TeV p-p collisions provided by the LHC. Clearly the 
different characteristic bands for various particles including deuterons and tritium 
are visible for particles and anti-particles, respectively. In order to achieve PID in a 
certain momentum region, histograms are filled and fitted with multiple Gaussians. 
An example is shown in Fig. [14] For a given ionization signal in this momentum 
region one can derive the particle yields, which represent the probabilities for the 
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Figure 14: Distribution of the difference between the measured ionization signals 
and the one expected for Kaons for a momentum slice of 50MeV/c width. The 
lines are fits, indicating that a sum of four Gaussians represents well the data. The 
peak centered at zero reflects the abundance of kaons, the other peaks represent 
other particle species. 



particle to be a kaon, proton, pion or electron. 

The ALICE TPC is in nearly continuous data taking mode since more than one 
year, providing tracking and PID information while the LHC delivers p-p collisions 
at ^ = 7 TeV. In November 2010 the first Pb-Pb collisions at ^/s]^ = 2.76 TeV 
were recorded. The energy resolution of the TPC in these heavy-ion collisions was 
found to degrade slightly to an average of about 5.3%. This effect is expected due to 
overlapping signals from neighboring tracks. The energy resolution for the highest 
multiplicity events falls off by 20% with respect to low multiplicities. Part of the 
increase is also explained by baseline fluctuations in the electronics due to large 
hit densities in single channels. These fluctuations can in the future be minimised 
using the signal tail cancellation and baseline correction features available in the 
TPC front-end electronics. In the course of the year 201 1, higher luminosities will 
be reached and space charge effects and gating efficiency will have to be accurately 
evaluated to avoid negative effects on the tracking and ionization sampling. 

3.8 Silicon detectors at LHC 

At the LHC, the signal amplitude information of the silicon strip (ALICE, AT- 
LAS and CMS) and drift (ALICE) detectors of the tracking systems make possible 
hadron ID via ionization measurements. However, since in solid media due to a 
stronger density effect the Fermi plateau is only slightly above the minimum, good 
tt/K (p/K) separation is achievable only up to 450MeV/c (up to 1 GeV/c). Again, 
the energy loss is estimated as a truncated mean in order to minimize the influence 
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of Landau fluctuations. In ALICE, the resolution of the measurement in two layers 
each of silicon strip and drift detectors is about 11%. In the CMS tracking system, 
the signal amplitudes from ten layers of silicon strip detectors may be used, but a 
resolution was not known to the author at the time of writing this review. 



3.9 Developments for future TPCs 

A TPC at a future linear collider will very likely be read out by micro-pattern 
devices like GEMs |fT9ll or MicroMegas ll34l . Their advantages are the high rate 
capability and a low ion feedback. Large TPC protopypes with different readout 
options were built tested [35]. The conventional readout of the micro-pattern de- 
vice, with mm-sized cathode pads, may actually be replaced by a digital device 
with pixel sizes of the order 55 /Lim. A family of ASIC chips is particularly suited 
for this purpose: MEDIPIX [36] and its successor TIMEPIX [37], which records 
also the time. This kind of setup has been shown to detect tracks of minimum 
ionizing charged particles (MIPs) with excellent single-electron efficiency and un- 
precedented spatial resolution. In this way, the cluster counting method may be 
used to assess the energy deposit through ionization: rather than measuring the de- 
posited charge, the number of primary collision clusters are counted. This avoids 
the problems related to the fluctuation of the energy transfer in single collisions 



(see Sects. [33] and |3.4| ) and should provide the ultimate resolution of ~2% ||38l . 
Concerning the implementation in a large detector like a TPC, the challenges will 
include equipping a large active area with the pixelized readout and handling huge 
amounts of readout channels. 



4 Time-of-flight 

Time-of-flight (TOE) measurements yield the velocity of a charged particle by 
measuring the particle flight time t over a given distance along the track trajectory 
L. The particle velocity fi = v /c = L/tc depends on its mass m and momentum 
p through 

/5 = . \ • (12) 
Thus, one can calculate the mass m from measurements of L, t and p: 



™ = f\/^-l. (13) 



4.1 Time resolution and separation power 

If two particles with masses uia and m^, respecively, carry the same momentum, 
their flight time difference can be calculated as 
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Figure 15: Particle separation with TOF measurements for three different system 
time resolutions (aroF = 60, 80 and lOOps) and for a track length L = 3.5m. 
Infinitely good precisions on momentum and track length measurements are as- 
sumed. 
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With p » mc the approximation ^/l + {mc/p)^ 
and with Eq. [4] the separation power becomes 



1 + (?Tic)^/2p^ can be used, 
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\tA - tsl 
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(15) 



aTOF '^p'^cttof' 
Here axoF is the resolution of the TOF measurement. Misidentification of par- 
ticles occurs at higher momenta, where the time difference \tA — tsl becomes 
comparable to (Ttof- Assuming a time resolution of 100 ps (60 ps) and requir- 
ing a separation of na-r^^p, = 3, the upper limits for the momentum are 2.1 GeV/c 
(2.7 GeV/c) for K/vr separation and 3.5 GeV/c (4.5 GeV/c) for K/p separation (see 



Fig. 15 1. A lower momentum threshold is defined by the curvature of the tracks 
in the magnetic field. Assuming a path length L = 3.5 m and a magnetic field of 
0.5 T, only particles with a momentum larger than about 300 MeV/c reach the TOF 
wall. 



4.2 Errors affecting the resolution 

The mass resolution of a TOF measurement is given by |[39]| : 



dm 

m 



dp 

P 
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dt dL\ 

t + t) 



(16) 
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Since in most cases 7 » 1, the mass resolution is predominantly affected by the 
accuracies in the time and length measurements, rather than by the accuracy of 
the momentum determination. In general the track length L and the momentum 
p can be measured with rather good accuracies in the tracking system {cTp/p ~ 
1%, ai/L 10^'^). The time-of-flight t is measured with a certain accuracy 
gtof by a dedicated detector system. A TOF measurement always consists of 
two time measurements: t = ti — to. The event start time to is measured with 
accuracy atg. The stop time ti > to is measured with accuracy at^ by the actual 
TOF detector, usually a large area detector system some distance away from the 
interaction region. The overall TOF resolution is then given by (t^qf ~ ^to + ^ti ■ 

4.3 Resistive Plate Chambers for time-of-flight measurements 

In order to optimise a TOF system the resolution axo F has to be improved. In cur- 
rent experiments for subnuclear research, detectors using scintillation (the property 
of luminescence when excited by ionizing radiation) are very commonly used due 
to their good time resolution. However, as the experiments have gradually in- 
creased in size, and since the readout of scintillator detectors is expensive, the cost 
per area became an important issue. The Resistive Plate Chamber (RPC) was found 
to provide a simpler technology at a much lower price. RPCs can be seen as the 
successor of the metalUc Parallel Plate Avalanche Counters, overcoming their dis- 
advantage of being very vulnerable to the destructive energy released by possible 
discharges or sparks. The electrode plates of RPCs aie made from a resistive ma- 
terial, which effectively restricts excessive signals or discharges to a well localised 
area of the detector. As a side effect, all charges cause a local drop of the elec- 
tric field in the gas gap and in this spot the detector becomes insensitive to further 
traversing particles for a time of the order of the relaxation time r = psoSr, where 
eo is the dielectric constant and p and £r are the volume resistivity and the relative 
permittivity of the resistive material. Even though the remaining counter area is 
still sensitive to particles, a limit to the overall rate capability is introduced through 
this effect. In general, accurate time measurements are feasible up to charged par- 
ticle rates of a few kHz/cm^, which is acceptable for many purposes (compare to 



Sect. 4.7 1. 



The RPC design was improved by introducing the concept of multiple gaps 
(Multi gap Resistive Plate Chamber, MRPC [40]), resulting in an increase of de- 
tection efficiency. Reducing the individual gap sizes HTI . while at the same time 
increasing the electric field, resulted in what is called the Timing RPC, a low cost 
detector with a time resolution comparable to that of scintillators. The space reso- 
lution can be chosen by different sizes of the readout strips. 

4.4 Detector physics of Resisitive Plate Chambers 

The physics of avalanches in RPCs is very complex, even though the simple ge- 
ometry of the device would not suggest so. During the development of large area 



32 



RPC systems for the LHC experiments, the detector physics of RPCs was there- 
fore studied in detail ||42|42lll|45l|46l|47l|48l|4l. The good detection efficiency 
of (single gap) RPCs is explained by the very large gas gain. A very strong space 
charge effect reduces the avalanche charges by many orders of magnitude. Detailed 
simulations show that the space charge field inside the avalanches actually reaches 
the same magnitude as the applied electric field. 

The very good time resolution of RPCs is due to the strong uniform electric 
field, which provokes the avalanche process immediately after primary ionization 
is deposited in the gas volume. As a consequence, the intrinsic detector time resolu- 
tion (no electronics) is determined by the avalanche statistics, and can be estimated 
using a simple formula (for one gas gap): 

1.28 

0"/ntr ~ 7 ^ • (17) 

[ttg - r}g)VD 

Here Ug and rjg are the Townsend and attachment coefficients of the gas mixture 
and fD is the the drift velocity of electrons. Typical operational parameters ll42]| 
for Timing RPCs are ag = 123 mm~^, -qg = lO.Smm^^ and vd = 210/imns~^. 
For these values an intrinsic time resolution of aintr = 50 ps is obtained for a 
single gap. By increasing the number of gaps n, the resolution can be improved, 
but for small values of n the improvement is not following a simple 1 / y/n scaling, 
because it is dominated by the gap which has the largest signal B2l . For n = 10 
gaps a resolution of ajntr = 20 ps is feasible. 

The drop of the electric field in the gas gap at high particle rates affects effi- 
ciency and time resolution. It may be calculated in a simple way by assuming the 
particle flux to be a DC current that causes a voltage drop when it passes through 
the resistive plate. For a single gap of width h the average field reduction is: 

(Ai?) = p^^Q , (18) 

where p and a are the volume resistivity and thickness of the resistive electrode 
material and ^> is the particle flux per area. Q is the avalanche charge, which of 
course depends itself on the field reduction, and has thus to be found in an iterative 
procedure. An approximate analytic expression for the field fluctuations around 
(A^) is given in Ref. 

4.5 ALICE time-of-flight detector 

In specific momentum ranges, the PID signals from the ALICE TPC can not be 
used to separate kaons, pions and protons, because the ionization curves are very 
close or overlap (see also Fig.[5]). The NA49 and STAR experiments have shown ll50l 
|5T]| how the hadron identification capability may be extended by combining the 
PID informations from TPC and TOF. Moreover, a reliable identification of elec- 
trons is possible up to a few GeV/c in this way, since the hadrons (which for a given 
momentum are slower) can be eliminated with a cut on the particle velocity. 



33 




Figure 16: Cross-section of an ALICE MRPC IHl. 



The ALICE TOE detector 1521 is based on MRPCs. 91 MRPCs each are assem- 
bled into gas tight modules in 1 8 azimuthal sectors (supermodules) surrounding the 
ALICE TPC and TRD (see Eig.[8]). The supermodules have about 9 m length, cov- 
ering the pseudo-rapidity interval \r]\ < 0.9. More than 150000 readout pads with 
size 2.5x3.5 cm^ cover an area of about 150 m^. The gas mixture in the ten 250 fim 
thin gas gaps is 90% C2H2 E4, 5% C4H10, 5% SEg. A cross-section through one 



MRPC is shown in Eig. 16 



4.5.1 Time resolution 

The design goal for the ALICE TOE was to reach a global TOE resolution of 
axoF ~ 100 ps. The following contributions have to be considered: 

C^TOF = 0-?o + ^tl = + ^"intr + (^Elec + ^^Clock + ^Cai • (19) 

Here at^ is the uncertainty in the absolute time of the collision, aintr is the in- 
trinsic time resolution of the detector technology (~20ps) and aEiec combines the 
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ALICE-TOF resolution measured with cosmic rays 
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Figure 17: TOF measurements with a cosmic data sample from the year 2009. 
Since the start and stop time were both measured with the ALICE TOF, the reso- 
lution of a single time measurement is obtained from the histogram by dividing its 
width by \f2. The resulting resolution of 88.5 ps corresponds to ut^ from Eq. 
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intrinsic time jitter of the amplification electronics (~20 ps) and of the time-to- 
digital conversion (~30ps). aciock is an uncertainty arising from the distribution 
of the digital clock from the LHC to the experiment (~15 ps) and through the elec- 
tronics chain (~10ps). Finally, a contribution acai includes all effects that can 
be parametrised and that are thus accessible to calibration methods in order to be 
minimised. These are mainly 

1. the relative timing of the different channels given by cable lengths {time 
offset), 

2. the influence of the finite rise time of the ampliflying electronic^(?/me slew- 
ing) and 

3. the influence of the impact position of the particle through signal propagation 
delays on the pad {time walk). 

It is forseen that the contribution acai can be kept below 30 ps f\E\ . Already with 
the limited statistics available from cosmics data samples before the start of the 
LHC it was possible to reach a value of acai = 70 ps. The resolution at^ with this 
"intermediate calibration" (compare to Tab. [2]) is shown in Fig. 17 



'in tlie ALICE TOF front-end electronics the channel-to-channel time slewing corrections are 
measured by a Time-Over-Threshold circuit, which provides a time width approximately propor- 
tional to the signal charge. 
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The absolute time of the colhsion to is fluctuating with respect to the nominal 
time signal from the LHC according to a Gaussian distribution with a cr of about 
140 ps. This is due to the finite length of the colliding bunches of particles. In 
ALICE, to is precisely measured by a dedicated detectoi[^(TO) |[T3^ . which reduces 
this uncertainty to at^ = 50 ps. 

Table 2: Overall TOP resolution together with the two main contributions cTty and 



O'TOF 


o-to 




Remark 


100 


50 


70 


proton collisions, intermediate calibration 


80 


50 


30 


proton collisions, optimal calibration 


65 


10 


30 


Pb-Pb collisions, optimal calibration 



With the much larger statistics that is available from data samples with col- 
lisions provided by the LHC, the work now focuses on reducing the contribution 
acai- A further improvement can be achieved for high multiplicity Pb-Pb events, 
where the accuracy of the event start time measurement is expected to improve to 
atg = 10 ps. In that case an overall system resolution axoF = 65 ps is possible 
(compare to Tab.|2]l. 



4.5.2 PID performance and outlook 



Fig. 18 shows the correlation between the track momentum measured by the ALICE 
tracking detectors and the velocity /3 as measured by the TOE system. The data 
were recorded during the first = 7 TeV p-p collisions in 2010 and clearly show 
two times four bands corresponding to four particle species and their corresponding 
anti-particles. Data outside these bands indicates errors in the association of tracks 
to TOE PID signals. A particle is identified when the corresponding point in the di- 



agram of Eig. 18 can be associated with only one theoretical curve (Eq. 12 1 within 
the measurement errors. Eor certain momentum regions histograms are filled and 



fitted with multiple Gaussians, similar to the procedure described in Sect. 3.7 An 



example is shown in Eig. 19 Eor a given TOE measurement one can derive the 
particle yields, which represent the probabilities for the particle to be a kaon, pro- 
ton, pion or electron. The power of combining TOE and ionization measurements 



for PID is demonstrated for in Eig. 20 



4.6 Other TOF detectors 

RPCs have found applications in many particle physics experiments in many vari- 
ations, as can be seen in the proceedings of the RPC workshop series |[53l l54l l55l 

'"Alternatively, the information from the TOF detector itself can be used to determine to in events 
with at least three tracks with an associated TOF signal by means of a combinatorial algorithm. 
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Figure 18: Velocity /? = w/c as measured with the ALICE TOF detector as a 
function of the particle momentum p multiplied with the particle charge number Z 
for a data sample taken with ^/s = 7 TeV collisions provided by the LHC in the 
year 2010. No data is available for momenta <300 MeV/c, since these particles do 
not reach the detector due to the curvature of their tracks in the magnetic field. 
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Figure 19: Histogram of the TOF signal for a given momentum window of 
lOOMeV/c width. The data sample was recorded with ^/s = 900 GeV p-p col- 
lisions provided by the LHC. The lines are fits to the data, indicating that a sum of 
three Gaussians well represents the data. 
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Figure 20: Demonstration of ALICE'S PID capabilities by combining the ioniza- 
tion measurements in the TPC and the mass calculated using the TOF signal. Elec- 
trons, pions, kaons and protons are clearly visible in a wide momentum range. The 
data sample was recorded with = 900 GeV p-p collisions provided by the 
LHC. 

I56l l57l l58l . The STAR experiment has installed a large area (64 m^) TOF sys- 
tem similar to the ALICE TOF in their barrel, surrounding the TPC (5911601. The 
modules have six gas gaps of 220 fim width and are filled with a gas mixture of 
95% C2H2F4 and 5% i-C4Hio. The STAR experiment also plans to install a Muon 
Telescope Detector (MTD) based on Timing RPCs to identify muons |[6TI . The 
FOPI experiment has upgraded their TOF system in 2007 by installing Multistrip 
Multigap Resistive Plate Chambers (MMRPCs) with 8 gaps of 220 /um |62l. For 
the HADES experiment a similar upgrade of the inner TOF wall is planned for the 
near future ||63]| . 

4.7 Developments for future TOF systems 

Driven by the environment in which TOF systems are planned to be operated in 
future experiments, a number of R&D projects currently address the rate limitation 
of Timing RPCs. The requirements for the TOF wall of the CBM experiment, 
for example, include a rate capability up to 25 kHz/cm^ with a time resolution of 
80 ps on a large system (120m^) |[64l . The effective drop of the operating voltage 
at larger rates (see Sect. |4.3| ) can be reduced either by lower electrode resistivity 
or by thinner electrode plates. Different electrode materials are being explored, 
including low-resistive phosphate and silicate glass [65, 66 1. With electrodes made 
from ceramics the practical feasibility of accurate timing measurements with RPCs 
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Figure 21: Schematic drawing of a TOF counter based on a MCP-PMT ITTOII . The 
MCP-PMT actually detects Cherenkov photons emitted in the quartz glass entrance 
window and has a very good time resolution. 

at rates up to 500 kHz/cm^ was established in test beams Il67ll68]| . 

TOF measurements with charged particles can also be carried out using Micro- 
Channel Plate (MCP) Photo Multiplier Tubes (PMTs, see Fig. |2T])), which allow 
time measurements with a few ps precision. The MCP-PMT actually detects Che- 
renkov photons emitted in the entrance window of the PMT (for details on Cheren- 
kov radiation see Sect. |5]l ||69ll . Time resolutions of at^ w 6ps (with an intrinsic 
resolution of the detector of ajntr < 5 ps) were observed in test beams ifTOHTTI . At 
this scale, the resolution of the measurement of the event start time becomes 
dominant and must be minimised as well. It remains to be shown that this tech- 
nology can be implemented on large surfaces at an affordable price and can thus 
compete with MRPCs. 

5 Cherenkov imaging 

In 1958 P. A. Cerenkov, I. Y. Tamm and I. M. Frank were awarded with the Nobel 
Prize in physics for the discovery and interpretation of the Cherenkov effect. Che- 
renkov radiation is a shock wave resulting from a charged particle moving through 
a material faster than the velocity of light in the material. The Cherenkov radiation 
propagates with a characteristic angle with respect to the particle track Gc, that 
depends on the particle velocity: 

cos(ec) = , (20) 
pn 

where n is the refractive index of the material. In general, the refractive index 
varies with the photon energy: n = n{E) (chromatic dispersion). Since | cos(©c)| < 
1, Cherenkov radiation is only emitted above a threshold velocity Pt = 1/"- and 

7* = 1/(1-/3?)^/^. 

In general, Cherenkov detectors contain two main elements: a radiator through 
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which charged particles pass (a transparent dielectric medium) and a photon de- 
tector. The number of photoelectrons (A'p.e.) detected in a given device can be 
approximated as ||23]| 

TVp.e. «A^o^'isin2(ec), (21) 

where L is the path length of the particles through the radiator, ze is the particle 
charge and A'^o is a quantity called the quality factor or figure of merit. As Che- 
renkov radiation is a weak source of photons, the light transmission, collection 
and detection must be as efficient as possible. These parameters are contained in 
Nq, as well as the photon collection and detection efficiencies of the photon detec- 
tor. Typical values of Nq are between 30 and 180cm^^. Three different types of 
Cherenkov counters can be distinguished: 

1 . Threshold counters measure the intensity of the Cherenkov radiation and are 
used to detect particles with velocities exceeding the threshold f3f A rough 
estimate of the particle's velocity above the threshold is given by the pulse 
height measured in the photon detector. 

2. Dijferential counters focus only Cherenkov photons with a certain emission 
angle onto the detector and in this way detect particles in a narrow interval 
of velocities. 

3. Imaging Cherenkov detectors make maximum use of the available informa- 
tion (Cherenkov angle and number of photons) and can be divided in two 
main categories: RICH (Ring Imaging CHerenkov) and DIRC (Detection of 
Internally Reflected Cherenkov light) devices. 

In this paper the focus lies on RICH devices, since the detectors used at the 
LHC experiments and discussed here all fall into this category. 

5.1 Cherenkov ring imaging 

RICH detectors resolve the ring shaped image of the focused Cherenkov radiation. 
From the knowledge of the particle momentum p and Cherenkov angle 0c a de- 
termination of the mass of the charged particle is possible. Combining Eqs.[T2|and 
|20]yields 

m = - \/ in? cos^iQc) — 1 • (22) 
c 

In a RICH device, the Cherenkov radiation is emitted in the radiator and collected 
by a photon detector, usually after being transmitted by optical means. The first 
RICH was developed by A. Roberts in 1960 11721 . The particular design was lim- 
ited by small angular acceptance and surface area and by low quantum and single 
electron counting efficiencies. These problems were overcome by J. Seguinot et 
al. in 1977 [i73il : A wire chamber was converted into an efficient single photon 
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detector by adding a photosensitive molecule into the gas mixture and replacing 
one of the cathode planes by a wire mesh and a UV transparent window. 



5.1.1 Radiator 

A RICH detector is designed to measure velocities in a specified momentum range 
by using a Cherenkov radiator with refractive index n chosen such that the Che- 
renkov angle varies with velocity, from threshold to the highest anticipated mo- 
mentum. The thickness L of the radiator is adjusted in order to assure a sufficient 



number of photoelectrons for the given momentum range (see Eq. 21 1. Two com- 



mon RICH designs can be distinguished, according to the value of n: 

1. If a dense medium (large n) is used, only a thin radiator layer (~1 cm) is 
required to emit a sufficient number of Cherenkov photons. The photon 
detector is located some distance away behind the radiator {expansion gap, 
usually about 5 to 10 cm), allowing the light cone to expand and form the 
characteristic ring-shaped image. A detector designed in such a way is called 
proximity-focusing (i.e., the focusing is achieved by limiting the emission 



region of the radiation). An example is shown in Fig. 22 



If a gaseous medium (n 1) is used, particles have to pass a thicker layer 
(>50cm) in order to emit a sufficient number of Cherenkov photons. In 
general, fiuorocarbon gases are chosen because they have a low chromatic 
dispersion (i.e. n does not depend strongly on the photon energy). The light 
is focused by spherical or parabolic mirrors onto the photon detectors where 



ring-shaped images are formed. An example is shown in Fig. 24 



5.1.2 Optics 

The quality of the optics of a RICH system influences the precision in the Che- 
renkov angle measurement and the figure of merit Nq. Mirrors should have high 
reflectivity to avoid photon loss. Traditional constructions use a glass substrate 
with coating of Al for the reflective surface and MgF2 or Si02 for protection. Such 
mirrors achieve values for the reflectivity of ~90%. For applications where a mini- 
mum of material budget is required, mirrors based on carbon fibre or Be substrates 
are used (see e.g. ||74]| and |[75]| ). 



5.1.3 Photon detection 

Cherenkov photons are converted to photoelectrons in photocathodes, typically 
made from Csl or bialkali with low work functions. The photon detection effi- 
ciency for a given photon detector is the product of the quantum efficiency (the 
probability that an incident photon produces a photoelectron, typically 20 to 30%) 
and the collection efficiency (the efficiency for detecting the photoelectron, typi- 
cally 80 to 90%). Photon detectors for RICH systems generally fall in two cate- 
gories: wire chambers with Csl coated cathode pads or detectors based on vacuum 
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tubes with bialkali photocathodes. Examples of both categories are described in the 
following (Sects. 5.3 and 5.4 1. Low noise is a general requirement to ensure the 
detection of single photoelectrons. Recently also the readout speed of the devices 
became increasingly important in order to cope with the highest event rates. 



5.1.4 Pattern recognition 

In the busy environment of hadronic collisions (such as at the LHC) many tracks 
may pass through the detector, which leads to overlapping rings. The assignment 
of photon hits to rings and the association of rings to tracks requires pattern recog- 
nition. Most approaches rely on the use of the particle track prolonged from the 
tracking system to seed the ring search. After transformation through the optics 
of the RICH, the track will lie inside the ring. For RICH designs where the track 
is known to always point through the ring center the ring search is actually quite 
simple and corresponds to the search for a peak in the number of photon hits versus 
radius from the track. In any case, the ring search efficiency needs to be studied in 
detail using Monte Carlo methods. 



5.2 Angular resolution and separation power 

The particle velocity /? can be calculated from the reconstructed Cherenkov angle 



according to Eq. 20 The resolution of the velocity measurement is given by 



^ = tan(ec) o-oc , (23) 

where cjoc is the resolution of the Cherenkov angle measurement. In practical 
counters uq^ varies between about 0. 1 and 5 mrad depending on the size, radiator 
material and length, and on the photon detector. The Cherenkov angle is deter- 
mined by A'p.e. measurements of the angles of emission of the single Cherenkov 
photons. With the average angular resolution for the angle measurement for the 
single photoelectron ue . , the total resolution becomes 



al^=\^ \ +4zo.. (24) 



The term aciob combines all contributions that are independent of the single pho- 
toelectron measurement. These include misalignment of the photon detector mod- 
ules, the resolution deterioration due to multiple scattering and background hits 
and errors in the calculation of the reconstructed track parameters. For the single 
photoelectron angular resolution, the following contributions have to be consid- 
ered: 

C^e, = 0-|p + 0"C/iro + C^Det J (25) 

where 
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1. (Tep is the geometrical error related to the emission point. For proximity- 
focusing devices the emission point is not known, and all photons are in 
general treated as if emitted at one point of the track through the radiator. 
This leads to some smearing of the reconstructed angle. Also for other de- 
signs the geometry of the optics (e.g. the tilting of the focusing mirror) can 
lead to a dependence of the image of a Cherenkov photon on its emission 
point on the track. 

2. achro is the error due to the chromatic dispersion of the radiators. Cherenkov 
photons are emitted with a distribution of wavelengths. Since the refractive 
index varies with wavelength, the Cherenkov angle is also spread. 

3. aDet is the error introduced by the finite granularity of the detector (e.g. pixel 
or pad size), which limits the precision. 

These contributions depend on the signal-to-noise ratio (detector gain) and can 
be a function of the impact angles of the track with respect to the detector surface. 

Assuming two particles of momentum p and mass niA and ms, respectively, 
then the Cherenkov angles @c,A and Qc,B are measured. The resolution 
determines the momentum range over which the two particles can be separated, 
since the Cherenkov angle saturates and thus the ring radii get closer. Using Eq. [4j 
the separation power can be calculated as 

For velocities /3 1 well above threshold , the separation power can be approx- 
imated as ll23l 

'm%-m\\. (27) 



Note the similarity to the formula for the separation power using the TOF tech- 
niqu e (Eq. [T5] ). In the case of a RICH there is however an additional factor of 
— 1, which allows to adjust the detector configuration in order to achieve 



the desired momentum coverage. Together with Eqs. 21 and 24 it is clear that the 
best separation is achievable by designing counters able to detect the maximum 
number of photons (large A'^o) with the best single photoelectron angular resolution 
(small as,)- 



5.3 ALICE HMPID detector 

The ALICE High-Momentum Particle IDentificaton (HMPID) 1771 was designed 
to enhance the PID capabilities of the ALICE experiment. The HMPID is a RICH 
system optimised to extend the useful range for K/tt and K/p discrimination up 
to 3 GeV/c and 5 GeV/c, respectively. Moreover, light nuclei (alpha particles, 
deuterons, tritium) can be identified. The HMPID is a single-arm array with a 



43 



charged particle 




Figure 22: Left image: schematic layout of an ALICE HMPID module, sliow- 
ing radiator, expansion gap and photon detector |[T3l . Right image: example 
ring as seen by the HMPID. The ring has a radius of the order 10 cm, but its 
shape is distorted due to the impact angle of the particle. 



reduced geometrical acceptance of |r/| < 0.6 in pseudo-rapidity and 1.2° < < 
58.8° in azimuth. The detector is based on proximity-focusing RICH counters and 
consists of an array of seven identical modules of about 1.5 x 1.5 m^ each. 



5.3.1 Radiator 

In each module, the radiator is a 1.5 cm thick layer of C6F14 liquid (perfluorohex- 
ane) with an index of refraction of n 1.3 at A = 175 nm. The material has a low 
chromaticity, keeping small the corresponding term in the angular resolution. The 
expansion gap has 8 cm thickness. 



5.3.2 Photon detector 

Cherenkov photons are detected by wire chambers with a thin layer (300 nm) of Csl 
deposited onto the cathode pad plane. Csl is a stable alkali halide crystal with high 
quantum efficiency of 25% at a wavelength of 175 nm 117811791 . The total size of 
the photosensitive area is 11 m^, making the HMPID the largest scale application 
of the Csl technique so far. The expansion gap and wire chambers are filled with 
CH4 at atmospheric pressure. The size of the cathod pads is 8 x 8.4 mm^. 



5.3.3 Performance and outlook 

The identification of charged particles with known momenta in the HMPID re- 
quires their tracks to be extrapolated from the central tracking devices (ITS, TPC 
and TRD) and associated with an ionization cluster in the HMPID. For the partic- 
ular setup of a proximity-focusing RICH, the shape of the ring is determined by 



the incident angle of the particle track (see Fig. 22 1, which makes the Cherenkov 
angle reconstruction quite complex. With a single photoelectron angular resolution 
(Te, ^12 mrad for perpendicular tracks fSOl, the resolution in the reconstruction of 
the Cherenkov angle is cre^ w 3.5 mrad (for perpendicular tracks and for /? w 1). 
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Figure 23: Dependence of the Cherenkov angle measured by the ALICE HMPID 
on the particle momentum. The lines are the theoretical curves calculated using 
Eq. [20] with the refractive index n = 1.3. No data is available for momenta below 
about 600MeV/c, because only particles with higher momenta reach the detector. 



Measured Cherenkov angles for different particles as a function of momentum are 



shown in Fig. 23 In the high multiplicity environment of Pb-Pb collisions the 
detection of Cherenkov rings is more challenging and the angular resolution is 
slightly worse. However, at the time of writing this article no quantitative results 
from the first heavy-ion data taking period (November 2010) were available to the 
author yet. 

In order to extend the PID range up to around 30 GeV/c, it was proposed to add 
a Very High Momentum PID (VHMPID) detector to ALICE [81 ]. The physics re- 
quirements have driven the choice towards a RICH detector using a C4F10 gaseous 
radiator (n = 1.0014) with length L = 80 cm in a focusing configuration using a 
spherical mirror. 

5.4 LHCb RICH 



The LHCb experiment was introduced in Sect. 2.3 Important for the underlying 
physics goals is a K/vr separation in the momentum range 2 < p <100GeV/c, 
which is achieved by two RICH detectors US HI [83. RICHl and RICH2 are both 
characterised by a design that places the photodetectors outside the acceptance of 
the LHCb spectrometer in order to limit the degradation of the resolution of the 
tracking systems due to interaction with the detector material. A set of spherical 
and flat mirrors projects the Cherenkov light onto the detector plane. The spherical 
mirrors are placed inside the spectrometer acceptance and special care was taken 
to minimise the fraction of radiation length while at the same time ensuring the 
mechanical integrity of the optics. The spherical mirror of RICHl is made from 
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Figure 24: Left image: schematic layout of the LHCb RICHl detector. The 
acceptance of ib250mrad is indicated. Right image: schematic view of an HPD 
of the LHCb RICH system. Both figures are taken from Ref. [ 18,1 . 



carbon fibre, with thickness corresponding to only about 0,01 Xq. The schematic 



layout of RICHl is shown in Fig. 24 



5.4.1 Radiators 

In both detectors gaseous radiators are used. In RICHl a second, denser radiator 
made from silica aerogej^ is installed in addition. Together the three radiators 
cover the large momentum range that is required. The parameters of the three 
radiators are summarised in Tab. [3] The Cherenkov angles for the three different 
materials and for different particles as a function of the momentum are shown in 



Fig. 25 



5.4.2 Photon detectors 

For the photon detectors the requirements are very demanding: they should be sen- 
sitive to single photons with high efficiency, cover a large area (about 4m^) with 
a granularity of about 2.5 x 2.5 mm^, and have readout fast enough to match the 
LHC bunch-crossing separation of 25 ns. LHCb chose Hybrid Photon Detectors 
(HPD) [86 1, which combine the photon sensitivity of vacuum PMTs with the ex- 
cellent spatial and energy resolutions of silicon sensors. The conversion of the pho- 
toelectrons takes place in the quartz window with bialkali photocathode and they 
are then accelerated through a potential difference of about 20 kV and detected by 
a silicon pixel chip. HPDs are sensitive to single photons in the wavelength range 



"silica aerogel provides refractive indices that range between about 1.0026 and 1.26 1841 [85 
closing the gap between gases (n ~ 1) and liquids/solids (n > 1.3). 
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Table 3: Some parameters of the LHCb RICH detectors. The measured single 
photoelectron angular resolutions liSTl are for the preliminary alignment available 
from the first data sample with p-p collisions at ^/s = 7 TeV. 





RICHl 


RICH2 




Silica aerogel | C4F10 


CF4 


Momentum range [GeV/c] 


<10 


10 < p < 60 


16<p<100 


Angular acceptance [mrad] 


vertical 


±25 to ±250 


±15 to ±100 


horizontal 


±25 to ±300 


±15 to ±120 


Radiator length [cm] 


5 


95 


180 


Refractive index n 


1.03 (1.037) 


1.0014 


1.0005 


Maximum Cherenkov angle [mrad] 


242 (268) 


53 


32 


Expected photon yield at /3 ~ 1 


6.7 


30.3 


21.9 


[mrad] 


expected 


2.6 


1.57 


0.67 


measured 


~7.5 


2.18 


0.91 



from 200 to 600 nm with a quantum efficiency of about 31%. About 65% of the 
detection plane is covered with read out pixels. 

The photoelectron trajectories are deformed by the magnetic field of the large 
LHCb dipole magnet. This effect can be corrected by comparing reproducible 
light-spot patterns, that are illuminated onto the photon detector plane, with mag- 
netic field on and off. 



5.4.3 Performance 

A typical (simulated) LHCb event in the RICHl detector is shown in Fig. [26] The 
expected single photoelectron angular resolutions are given in Tab. [3] After the 
first alignment carried out with data from p-p collisions, the achieved resolutions 
are already approaching those values ||87]| . For the aerogel radiator, the refractive 
index was found to change from its initial value of n = 1.03 to n = 1.037 in the 
experiment, which is probably explained by absorption of C4F10, since the aerogel 
material is not sealed from the gas volume. As a consequence, the measured photon 
yield and angular resolution are slightly worse than expected. Using the expected 
values, the separation power ria-Q^ can be calculated for all three radiator materials 



according to Eq. 26 The results are shown in Fig. 27 

Particle identification with the LHCb RICH system uses a global pattern recog- 
nition and maximum likelihood method. All found tracks in the event and all three 
radiators are considered simultaneously and the observed pattern of hit pixels in 
the RICH photodetectors is matched to that expected from the reconstructed tracks 
under a given set of particle hypotheses. The likelihood is maximised by varying 
the particle hypotheses of each track to be an electron, muon, pion, kaon or proton. 
In this way efficient K/vr separation is achieved with high purity [ 1 8 1 : the average 
efficiency for kaon ID for momenta between 2 and lOOGeV/c is at the level of 
90%, with a corresponding average pion misidentification rate below 5%. 
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Figure 25: Cherenkov angles as a function of momentum for different particle 
species and for the three different values of the refractive index n corresponding to 
the three radiator materials used in the LHCb RICH setup. 



(cm) 
60 




Figure 26: Typical, simulated LHCb event in the RICHl detector lITSl . The data 
from the two photodetector planes are drawn in the upper and lower halves. 
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momentum (GeV/c) 

Figure 27: Particle separation achievable witli Cherenkov angle measurements for 
the three different radiator materials used in the LHCb RICH detectors. As Che- 
renkov angular resolutions the expected values from Tab. [3] were used. 

5.5 Other Cherenkov detectors 

Cherenkov counters have found applications in many particle physics experiments 
in many different configurations, as is documented in the proceedings of the RICH 
workshop series [88, 89, 90,, ,911 |92l |93l |94J. To demonstrate this variety, here a 
selection of other Cherenkov detectors ouside the LHC is given. 

• A RICH with a 16 m long He (n = 0.000035) radiator, 8 m focal length 
mirrors and wire chambers with triethylamine (TEA) as photosensitive com- 
ponent was used in a spectrometer at FNAL ll95l l96l and provided K/vr sep- 
aration up to 200GeV/c. 

• The HERA-B collaboration showed that a RICH with gaseous (C4F10) ra- 
diator and multianode PMT readout can be safely operated at high track 
densities, and no degradation of performance was observed in five years of 
operation [97 1. 

• The central region of the RICH of the COMR^SS experiment was upgraded 
by replacing the wire chamber based photon detector with multianode PMTs 
very similar to the ones employed in the HERA-B experiment ||98l . 

• In the Belle spectrometer a Cherenkov counter based on an aerogel radiator 
and fine mesh PMT readout is used in a threshold configuration ||99ll : The 
refractive index is chosen such that pions emit Cherenkov radiation while 
kaons stay below threshold. 

• At the CLEO-III experiment a proximity-focusing RICH counter based on 
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LiF as solid radiator and wire chambers with TEA as photosensitive compo- 
nent is installed IIIOOI . 

• The RICH of the HADES experiment employs a gaseous (C4F10) radiator 
and a gaseous photon detector with Csl converter and acts as a hadron blind 
trigger device lilOlL 

• The Hadron Blind Detector (HBD) of the PHENIX experiment is built for 
electron ID. It uses a gaseous radiator (CF4) and as photon detector three lay- 
ers of Gas Electron Multipliers (GEM) with Csl coating on the first layer [ 102J 

• The DIRC of the BaBar experiment [103| uses 4.9 m long, rectangular bars 
made from synthetic fused silica (average refractive index n 1.473) as 
Cherenkov radiator and light guide. The photons are focussed by a "pin- 
hole" and the image is expanded through a standoff region filled with 6000 1 
of purified water onto an array of PMTs placed at a distance of about 1.2 m 
from the bar end. 



5.6 Developments for future RICH detectors 

The Cherenkov counters of the BaBar and Belle spectrometers have contributed 
significantly to the recent advances in flavor physics. In both cases, upgrades are 
planned to further improve the performance and to allow them to work at higher 
event and background rates. The LHCb collaboration is planning to upgrade the 
detector around the year 2016, following first data taking phase II104II . The upgrade 
strategy involves increasing the design luminosity by a factor of 10 to about 2 x 
10^"^ cm^^s"^. In that case the full experiment is read out at 40 MHz (currently 
1 MHz). As a consequence the HPDs need to be replaced with a faster technology. 
There are a few promising candidates, however, here the developments in the wide 
area of photon detectors are not discussed. A summary is given for example in 
Refs. ifTMlfToa . 



5.6.1 Proximity- focusing aerogel RICH 

The drawback of the proximity-focusing RICH technology is that for thicker ra- 
diators (needed to increase the number of photons) the emission point uncertainty 



(aEP in Eq. 25 1 leads to a degradation of the single photon resolution. This lim- 



itation can be overcome by using a radiator made from several layers of aerogel 



with refractive indices gradually increasing, as is shown in Fig. 28 The corre- 
sponding Cherenkov rings are overlapping on the photon detector, which leads to 
an optimised angular resolution. Such a proximity-focusing RICH with inhomo- 
geneous aerogel radiator is proposed for the endcap region of the upgraded Belle 
spectrometer |[T08l[T09il . 
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Focusing type 




Figure 28: Schematic view of a proximity-focusing RICH with an inhomogeneous 
aerogel radiator in the focusing configuration. Cherenkov photons from the differ- 
ent radiator layers overlap, thus minimising the error in the angular resolution due 
to the uncertainly in the emission point H 10711 . 




charged particle 

Figure 29: Schematic side view of a TOP counter UllOi Cherenkov photons are 
guided to the photon detector by total internal reflections. The difference in the 
propagation time for two particle types (here kaons and pions) is due to the emis- 
sion angle of the Cherenkov radiation and can be used to enhance PID information. 



5.6.2 Time-of-propagation counters 

For the barrel region of the upgraded Belle spectrometer, a Time-of-propagation 
(TOP) counter is being studied HI 101 llllll . Here the two dimensional information 
of a Cherenkov ring image is represented by the time-of-arrival and impact po- 



sition of the photons at the quartz bar exit window of a DIRC (see Fig 29 1. The 
construction is more compact than a DIRC, since the large expansion volume is not 
needed. A similar concept combining TOF and RICH techniques, named TORCH, 
was proposed for the LHCb upgrade as a possible solution for the identification of 
low momentum hadrons (p < lOGeV/c) III 121 . It uses a large quartz plate to pro- 
duce Cherenkov light and identifies the particles by measuring the photon arrival 
times. 
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Figure 30: Schematic view of the AMS 02 spectrometer II113II . 



6 PID in space 

Non-accelerator experiments have become increasingly important in particle phy- 
sics. An impressive example of a space-based experiment using a number of fa- 
miliar PID techniques is the Alpha Magnetic Spectrometer (AMS 02 II113I . see 



Fig. 30 1. The purposes of the experiment are to search for cosmic antimatter and 
dark matter and to study the composition and energy spectrum of the primary cos- 
mic rays. The AMS 02 is set for installation on the International Space Station 
(ISS). It will be transported in the cargo bay of the space shuttle on its final mis- 
sion, scheduled for 2011. The experiment will be attached to the outside of the 
space station and gather data that will be sent to a ground station for analysis. 

The AMS 02 experiment is a state-of-the-art particle physics detector in space. 
The initial design included a superconducting magnet (0.9 T), which was however 
recently replaced by a permanent magnet (0. 14 T) which allows the experiment to 
remain operational for the entire lifetime of the ISS (up to 18 years). The setup 
further includes a high-precision, double sided silicon strip detector and 4 dedi- 
cated PID devices: a TRD for electron ID (straw tubes filled with a mixture of Xe 
and CO2 and fleece radiators), a proximity-focusing RICH with an aerogel radiator 
(n = 1.03) for determination of the particle velocities and charge^^ a TOF de- 
tector (four planes of plastic scintillators) providing velocity determination and a 



^In the RICH the charge of the incident particle is found by measuring the number of photoelec- 



trons using Eq. 21 
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fast trigger, and a sampling EM calorimeter (lead absorber plates and scintillating 
fibres). The whole system provides identification of particles and isotopes up to 
iron. 

7 Summary and outlook 

Particle identification (PID) is an important ingredient to particle physics experi- 
ments. Short-lived particles are reconstructed from their decay vertex and/or decay 
products. Some long-lived particles (leptons and photons) can be identified from 
the signatures they leave in the different layers of a typical experiment. Distin- 
guishing the different long-lived charged hadrons (pions, kaons and protons) is 
more challenging, but often their identification is crucial, e.g. for precision mea- 
surements of rare processes with high statistics. Usually dedicated detectors are 
required for the task, which are based on the determination of the particle's mass 
by simultaneous measurements of momentum p and velocity /3 = v/c. The ve- 
locity is obtained by one of four techniques: measurement of the energy deposited 
through ionization, measurement of the time-of-flight, imaging of Cherenkov radi- 
ation or detection of transition radiation (transition radiation is not covered by this 
review). An illustrative comparison of the methods is shown in Fig. |5] The Che- 
renkov imaging method is the most flexible, as it allows to tune the response of the 
detector by varying the refractive index (and the length of the radiator), and makes 
accessible also very high momenta (>100GeV/c). The Cherenkov and time-of- 
flight techniques can be used to close the gap in momentum in the MIP region, 
where PID by ionization measurements is not possible. All techniques are very 
well understood, with detectors being of the second, third or fourth generation, but 
the use and sophistication is still increasing as new generation detectors are built. 

At the LHC, PID is one of the main challenges. Different concepts are fol- 
lowed: the two specialised experiments, ALICE for heavy-ion collisions and LHCb 
for B physics, make use of specific techniques adapted to their physics program. In 
both cases PID at low transverse momentum (but high total momentum in the case 
of LHCb) is important. ALICE actually uses all existing PID techniques. The gen- 
eral purpose detectors on the other hand, ATLAS and CMS, follow the traditional 
setup of a particle physics experiment: a tracking system, EM and hadron calorime- 
ters and a muon system arranged in layers, with large acceptance for the identifi- 
cation of charged hadrons, photons and leptons. For their physics program, the 
challenges for precision calorimetry and precision muon detection are extremely 
high. Here the two experiments are actually following complementary approaches, 
but it can be expected that they are very competitive. All four experiments are very 
well advanced with the detector calibration. They have started to record data as the 
LHC began delivering p-p collisions in the year 2009. In the course of the year 
2010, the increasing luminosity allowed real stress tests of the systems. At the end 
of that year also the first heavy-ion collisions were recorded by ALICE, ATLAS 
and CMS. 
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